Signature of valley polarization in fractional flux periodicity of a graphene ring 
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We have studied the interplay of valley polarization and the Coulomb interaction on the energy 
spectrum, persistent current, and optical absorption of a graphene quantum ring. We show that the 
interaction has a dramatic effect on the nature of the ground state as a function of the magnetic flux, 
and that the absence of the exchange interaction for electrons in opposite valleys means that the 
singlet-triplet degeneracy is not lifted for certain states. The additional level crossings (fractional 
flux periodicity) due to the interaction directly leads to extra steps in the persistent current and 
intricate structures in the absorption spectrum that should be experimentally observable. By varying 
the width of the ring, the nature of the ground state at zero field can be varied and this is manifested 
in the measurable properties we discuss. 
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Graphene, a single layer of carbon atoms arranged in 
a honeycomb lattice, with its unusual electronic proper- 
ties, has claimed the center stage of condensed matter 
research for the past three years [l| . Theoretically inves- 
tigated sixty years ago 0] , interest in this system widened 
after free-standing graphene flakes were obtained exper- 
imentally by Geim et al. in 2004. The linear band struc- 
ture predicted for this material has been verified and sev- 
eral striking experimental observations have been made, 
including the 'half-integer' quantum Hall effect and a 
solid-state manifestation of the Klein paradox for mass- 
less Dirac fermions. Q 

In bulk graphene, the Fermi energy is located at the 
two inequivalent K points and the corresponding pairs 
of single-particle eigenstates in each of these 'valleys' are 
degenerate with each other. An experimental observa- 
tion Q of lifting of this valley degeneracy in high-field 
quantum Hall effect measurements has prompted several 
theoretical studies seeking the origin of this valley polar- 
ization The intriguing possibility of controlling the 
energy difference between electron states in opposite val- 
leys which would facilitate the idea of valley-tronics (uti- 
lizing the valley quantum number to control the system) 
would be an exciting development. 

Quantum rings of nanoscale dimensions are known to 
carry a persistent current: an equilibrium current driven 
by the magnetic field threading the ring. This is a direct 
consequence of the Aharonov-Bohm (AB) effect which 
manifests itself as periodic oscillations in the energy spec- 
trum of the electronic system as a function of the number 
of flux quanta entering the ring [|| . Impressive progress 
in fabricating nanosize quantum rings containing only a 
few electrons has led to equally notable results on the ob- 
servation of the properties of energy spectra, first noticed 
in Ref . @ , via magneto-absorption spectroscopy Q and 
in magneto-transport measurements @. The important 
role of the electron-electron interaction in this system was 
found to lift the degeneracy between states with different 
spin as a means to gain the exchange energy A di- 



rect consequence of this is the fractional flux periodicity 
that was indeed observed in subsequent experiments on 
semiconductor quantum rings containing only about four 
electrons (loj . 

Graphene rings have recently been fabricated and AB 
oscillations were observed in their conductance (ll| . The 
combined effect of the ring confinement and applied mag- 
netic flux is suggested theoretically to lift the orbital de- 
generacy arising from the two valleys in a controllable 
way [12j. Further, a ring with quantum point contacts 
has been shown to polarize the transport current with 
respect to the valley The important effects of the 

Coulomb interaction on the valley degeneracy and ground 
state properties have not yet been investigated however. 

In this work, we report on the effect of the Coulomb in- 
teraction on the energy spectrum, persistent current and 
optical absorption spectrum of a graphene quantum ring. 
We show that the interaction, the total valley quantum 
number and spin will dramatically change the nature of 
the ground state of a few-electron system. We find that 
the interaction causes drastic changes in the nature of the 
ground state as the flux varies and that the absence of 
the exchange interaction for electrons in opposite valleys 
means that the singlet-triplet degeneracy is not lifted for 
some states. The extra crossings in the spectrum which 
arc generated by the interaction manifest themselves as 
steps appearing in the persistent current (the fractional 
AB effect) and results in intricate structures in the ab- 
sorption spectrum. These effects are all experimentally 
measurable. Details of the periods of the oscillations of 
the persistent current depend on the width of the ring 
and hence the interplay between kinetic and Coulomb 
energies. 

We use the valley-symmetric from of the graphene 
Hamiltonian [12j 
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where V(r) is a mass term which describes the confine- 
ment of the electron, H. = v F (p ■ a) is the bulk graphene 
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Hamiltonian, a x z and t x z are Pauli matrices in the 

sublattice and valley spaces respectively, p = —iftX7 + eA 
and v F is the Fermi velocity. The vector potential is 
taken as A = ($>/2Tvr)e v> where $ is the total magnetic 
flux threading the ring. The index N stands for the pair 
of indices [to, t], where we assume the electrons are in the 
lowest Landau level, to is the orbital angular momentum, 
t = +1 in the K valley and r = — 1 in the K' valley. We 
write the wave functions for V{r) = using the dimen- 
sionlcss length p = r/W and energy e N — E N W/hv F 
where W is the width of the ring [Fig. [2^c)] as 
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-§-, b N is the normalization 



with p± = ^ ± |, to = r 
factor and sgn(x) = 1 for x > and sgn(x) = — 1 for x < 
0. The functions Hv (x) and H^\x) are respectively 
Hankcl functions of the first and second kind. 

The ring confinement is defined by the mass term [writ- 
ten as the potential V(r)] in the Hamiltonian. We em- 



ploy infinite mass boundary conditions [12], [14[ so that 
V(r) — ► oo outside the ring. This yields the boundary 
condition ip{p±) = r(n ■ o)ijj(p±). The coefficient a N 
is found by applying this condition at the inside edge, 
and b N is calculated numerically for each state and value 
of the flux via the normalization condition for the wave 
function. 

The interacting few-electron system is studied by 
adding the term 
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to the Hamiltonian where e g is the dielectric constant 
of monolayer graphene. The simple matrix structure of 
this operator is unchanged by the transformation to the 
valley-symmetric form. We evaluate (numerically) the 
matrix elements of this operator over the single-particle 
wave functions in Eq. ©. Using these single-particle 
matrix elements, we construct the many body Hamilto- 
nian and carry out an exact diagonalization procedure 
to determine the energy and eigenstates of the interact- 
ing system. The persistent current j is then calculated 
by taking the derivative of the ground state energy E 
of the few electron system with respect to the flux as 



To describe the absorption of incident light by the 
graphene ring we add a term to the Hamiltonian which 
describes the coupling of electrons to the field via the 
vector potential ^4 EM = 2A ctcos(k ■ r — tut). We assume 
that the radiation propagates as a plane wave with wave 
vector k, frequency u, and polarization described by the 
unit vector a. Then, the Hamiltonian can be written 



H = v F a ■ [p + eA 
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in the valley symmetric representation. The transition 
rate from state N to state TV' is calculated from 
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in the dipolc approximation. The integral (where we drop 
the coordinate dependence of the spatial functions for 
brevity) must be evaluated numerically. The intensity 
of the absorption is proportional to this transition rate 
and the area of the dots in the lowest panels of Fig.[lja), 
Fig. [TJb) and Fig. fjjc) scale with this quantity. In all 
figures we show the absorption of unpolarized light [i. e. 
a = (e x + e y )/y/2]. Equation ([6]) shows that transitions 
which change the angular momentum quantum number 
by ±1 are permitted, so long as the valley index remains 
the same. Where the initial or final states of the transi- 
tion are degenerate, we take the average of the intensity 
of all possible pairs of initial and final states. 

In Fig. [lja) we show the energy spectrum, persistent 
current and optical absorption for a single electron in the 
graphene ring with R/W = 10. The lifting of the valley 
degeneracy previously described causes the step in the 



persistent current at 



0.5. For < 6 < 0.5 



the ground state consists of one electron in the m = — i , 
r = — 1 state whereas for 0.5 < 4> < 1 the valley index is 
t = +1. For 0^0, transitions to the lowest-lying states 
to = +i, r = +1 and to = — |, r = —1 arc not allowed 
since the optical absorption cannot mix valleys. 

For two non-interacting electrons, the ground state 
consists of a pair of electrons with anti-parallel spins oc- 
cupying the same single-particle states as in the single- 
electron system [Fig. QJb)]. The persistent current re- 
flects the similarity between the ground states of the 
single-particle and N = 2 non-interacting system, and 
since there are now two electrons, the persistent current 
is doubled. The excited states can have varying degrees of 
degeneracy: If the quantum number pairs P = [to p ,t p ] 
and Q = [toq, Tq] of the two electrons are identical then 
there is only one permitted configuration of the electron 
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(a) Single electron (b) N = 2, no interaction (c) N — 2, Coulomb 
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FIG. 1: Energy spectrum, persistent current (middle pane) and optical absorption of unpolarized light (lower pane) by (a) a 
single electron, (b) two non-interacting electrons, and (c) two electrons with the Coulomb interaction included. States in the 
two-electron plots are labelled by the pair of quantum numbers (M, T) where M is the total angular momentum, and T is the 
total valley index. The area of the points in the absorption plots represent the intensity of the peak in arbitrary units. In all 
three plots, W = 10 nm and $ = 10, and e 9 = 2.5 [H]. 



spins, the singlet state. However, if P ^ Q then there are 
four degenerate possibilities: the singlet and three triplet 
states. 

When the Coulomb interaction is included [Fig. [He)], 
the picture changes drastically. To describe the two par- 
ticle states, we introduce the notation M = m 1 + m 2 for 
the total angular momentum and T = t 1 + r 2 for the 
total valley quantum number. The exchange interaction 
will split the degenerate singlet-triplet states when both 
of the electrons are in the same valley i.e. for T = ±2. 
In this case, the energy of the singlet does not contain 
any contribution from exchange and consequently has a 
rather higher energy than the corresponding triplet. This 
is exemplified by the (M = 0, T = 2) state. The triplet 
part experiences exchange and this reduces its energy 
sufficiently for it to form the ground state for <f> w 0.3 
with £ ~ 3.191. At the same flux the singlet state has 
e ~ 3.205 and is therefore not present in Fig.[TJc). On the 
other hand, the singlet and triplet parts of the (—1,0) de- 
generate state are not split by the exchange interaction. 

Because the Coulomb matrix elements depend on the 
angular momenta of the single-electron states involved, 
the size of the energy change will vary between different 
states. This is shown in the vicinity of the crossing of the 
lowest three states for <j> = 0. The M = ±1 states each 
contain two electrons in the same angular momentum 
state (m = ±i) so their interaction is stronger than the 
electrons in the M = state which has electrons in dif- 
ferent angular momentum states. This causes the ground 
state to become a degenerate singlet-triplet combination 
for a small range of flux. 

This intricate interplay of different-sized contributions 
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FIG. 2: The effect of the ring width on the ground state en- 
ergy, (a) The ■& = 3 curve (crosses) is plotted relative to 
the right-hand axis, and the ^ = 25 curve (circles) relative 
to the left-hand axis, (b) The energy difference between low- 
est singlet (e s ) and degenerate singlet-triplet (e D ) states at 
$/<I > o = 0. (c) The geometry of the ring. 



from the Coulomb interaction adds significant complex- 
ity to the ground state of the interacting system. In 
the non-interacting case the ground state is always com- 
prised of a singlet, but the interaction introduces several 
additional level crossings which give rise to ranges of the 
flux where the ground state becomes a triplet or degener- 
ate singlet-triplet state. Moreover, since the energy due 
to the Coulomb interaction depends on the angular mo- 
mentum of the state, the size of the ring will also be an 
important factor. In Fig. [2] (a) we plot the ground state 
energy for ^ = 3 (a wide ring) and -p^ = 25 (a narrow 
ring) to illustrate this dependence. The relative depth 
of the minima of the ground state energy vary, and the 
nature of the ground state at zero field changes from the 
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N = 3, no interaction N = 3, Coulomb 
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FIG. 3: (a) Non-interacting, and (b) interacting three electron 
energy spectrum, persistent current and optical absorption for 
-jS- = 10. Dashed lines denote two-fold degeneracy, solid lines 
four-fold degeneracy and thick solid line eight-fold degeneracy 
of the state. 

(—1, —2) singlet in a wide ring to the (0, 0) singlet-triplet 
in a narrow ring. This transition is revealed by the ab- 
sorption spectrum since the crossover to the degenerate 
ground state changes the spectrum to two closely-spaced 
low intensity peaks. In Fig. H^b), the difference in energy 
between the (0,0) and (—1,-2) is plotted as a function 
of -pp for a ring with R = 100 nm. The crossover for the 
ground state occurs at approximately -p^ = 7, indepen- 
dent of the value of R. 

For three non-interacting electrons in the ring, the 
ground state is composed of spin and valley unpolarized 
states (i.e. T = ±1). When the interaction is added, the 
contribution from exchange is largest for T — ±3 states 
so the low energy spectrum becomes much more compact, 
just as in the N = 2 case. Qualitatively, the effect of the 
interaction is the same as previously, so that the chang- 
ing nature of the ground state again demonstrates the 
complexity due to the absence of the valley degeneracy. 
However, because there are more possible combinations 
of states, the persistent current and absorption spectrum 
are correspondingly more complex in their structure. In 
particular it is not possible to have T — so the ex- 
change energy is always finite. However, its contribution 
is larger for T = ±3 states than for T = ±1 states. It is 
also the case that the width of the ring (and hence the 
relative strength of the interaction) will affect the detail 
of the ground state. 

To summarize, we have studied the effect of the 



electron-electron interaction on measurable quantities in 
a graphenc quantum ring. We find that the interplay 
of the interaction and the total valley quantum num- 
ber allow for an intricate manifestation of the breaking 
of valley degeneracy in this geometry. The change of 
the interacting ground state between singlet, triplet and 
degenerate singlet-triplet natures reveals the sensitivity 
of the exchange contribution to the total valley index. 
These changes in the ground state are manifested in the 
fractional nature of the AB oscillations in the persistent 
current, and in the steps and intensity changes in the 
absorption spectrum as the flux is varied. 
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